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Pierre Charland, a new DUG Member from Canada sent a couple of “Graphics bundles” 
containing about 200 DERIVE plots. This is a tiny mouth watering selection: 





What a wonderful design for a water glass. (It reminds me on Youth Style, Josef) 
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Dear DUG Members, 

this is not the first issue of 2007 but should be 
considered as last issue of 2006. Anyhow, it 
does not contain so many contributions but it 
shows an amount of 46 full pages. I didn't 
want to split Bjorn Felsager's article into parts 
and I wanted to have a series of three articles 
on a fine use of CAS in classroom. As you can 
see there are only two of them - Traffic Den- 
sity and Morley Triangle -, the third one about 
a pretty property of quartics must be postponed 
for the next DNL. When I was ready copying 
and pasting Peter Liike-RosendahTs paper 
Morley in the Mirror I had in mind to leave it 
for the next issue because the newsletter was 
full, but then I decided to keep it as a penitence 
for my delay delivering DNL#64 - and, to be 
honest I found it too nice to leave it unread two 
months in my computer. 

Johann Wiesenbauer prepared a new Titbits- 
contribution about semigroups and groups, 
(which will appear in DNL#65) and he prom- 
ised to revise his Titbits from earlier issues. 
Next revised version will be DNL# 13 with the 
first Titbits-column. (By the way, do you re- 
member Titbits# 1 with the Extended Euklidean 
Algorithm. 


There is no information in this DNL on the 
first page but instead of it you can find a first 
selection of DERIVE graphs which reached me 
some days before Christmas. Pierre Charland 
from Canada composed many many terrific 
graphs, some of them in black and white, some 
of them with only a few colours and some of 
them with plenty of colours. I'll insert a selec- 
tion of them in the next DNLs. 

Koen Stulens from Belgium announced a cou- 
ple of short articles on the use of TI-Nspire to 
demonstrate the features of this new CAS-tool. 

Finally I don't want to forget wishing you a 
healthy and successful year 2007. We are all 
looking forward what 2007 will unhide for the 
CASers from all over the world. Let's make 
the best of it. 

Best regards until March (or April?) 




Don't forget to order your copy of the Proceedings of DES-TIME 2004 in Dresden. 
It a CD with a contents of much more than 500 MB. 

Order at http : //shop . bk-teachware . com 


Download all DNL-DERIVE- and Tl-files from 

http : //www . austromath . at/ dug/ 
http : / / www . derive-europe . com/ support . asp?dug 
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The DERIVE-NEWSLETTER is the Bulletin of 
the DERIVE & CAS-77 User Group. It is pub- 
lished at least four times a year with contents 
of 44 pages minimum. The goals of the DNL 
are to enable the exchange of experiences 
made with DERIVE and the 77- 
%9l92ITitaniumlVoyage 200 as well as to cre- 
ate a group to discuss the possibilities of new 
methodical and didactical manners in teaching 
mathematics. 

As many of the DERIVE Users are also using 
the CAS- 77s the DNI tries to combine the 
applications of these modem technologies. 


Editor: Mag. Josef Bohm 
A-3042 Wurmla 
DXust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 

e-mail : noj o .boehm@pgv. at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to write 
their contributions in English to reinforce the 
international touch of the DNI. It must be said, 
though, that non-English articles will be 
warmly welcomed nonetheless. Your contribu- 
tions will be edited but not assessed. By sub- 
mitting articles the author gives his consent for 
reprinting it in the DNI. The more contribu- 
tions you will send, the more lively and richer 
in contents the DERIVE & CAS-77 Newsletter 
will be. 


Next issue: March 2007 

Deadline 15 Febmary 2007 


Preview: Contributions waiting to be published 

Two Stage Least Squares, M. R. Phillips, USA 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Another Task for End Examination, J. Lechner, AUT 

Tools for 3D-Problems, P. Lilke-Rosendahl, GER 

ANOVA with DERIVE & 7/, M. R. Phillips, USA 

Financial Mathematics 4, M. R. Phillips 

Hill-Encription, J. Bohm 

Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE, J. Bohm 
Henon & Co, J. Bohm 

Are all Bodies falling equally fast, J. Lechner 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Lilke-Rosendahl, GER 

Mathematics and Design, H. Weller, GER 

Diophantine Polynomials, D. E. McDougall, Canada 

Challenger Matrix Problems, G P Speck, New Zealand 

MuPAD, a special CAS-Guest, Th. Himmelbauer, AUT 

Graphics World, Change Count, P. Charland, CAN 

Precise Recurring Decimal Notation, P. Schofield, UK 

Problems solved using the Tl-Nspire, K. Stulens, BEL 

and Setif, FRA; Vermeylen, BEL; Leinbach, USA; Koller, AUT; Baumann, GER; 
Keunecke, GER, and others 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, DXust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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From: DERIVE computer algebra system 

On Behalf Of Fred J. Ty deman 

Sent: Thursday, December 28, 2006 8:31 AM 

To : DERIVE-NEWS0 JISCMAIL . AC . UK 

Subject: 7.5 digit log 

How do I find the 7 digit number in the range of 10 to 20 that has a base- 
10 log that is closest to the midpoint between two 7 digit numbers? 

For example, 

log(10. 00706, 10) is 1 . 00030 6_5037 .. . 
log (10 . 00918, 10) is 1 . 000398_4 994 . . . 

Fred J. Tydeman Tydeman Consulting 


In your example the midpoint is exact. 
10.00812 
Jim FitzSimons 


December 29 


I believe not: 

#1: LOG (10 . 00812, 10) 

# 2 : 1 . 000352 50402203308649343332674562 


December 29 


Fred J. Tydeman 

December 29 

10 A ( (LOG (10 .0070 6, 10) +LOG (10 . 00 918 , 10 )) /2 ) = 

SQRT (250406162027 ) /50000 is the midpoint on the log scale. 

This is 10.00812 rounded to 7 digits. 

Jim FitzSimons 


While that may be, it is not what I am asking for. 


December 29 


I am asking for an exact 7 digit number (not a rounded number) whose log is 
closest to the midpoint between two 7 digit numbers . 


Fred J. Tydeman 


December 29 

10.00812 is the exact midpoint. That is a fact. 

The logarithm function is one to one so the midpoint on the linear scale is 
in the middle on the log scale, but not exactly the midpoint on the log 
scale . 


You have to explain the problem better. 

I do not understand what you are asking. 


Is 10188730000 what you are looking for? 

That is not between 10 and 20, but that is not possible given 7 digit num- 
bers 10.00706 and 10.00918. 


[LOG (1018872 0000, 10) , 10 . 00812 , LOG ( 10188730000 , 10 ) ] 
[10.00811962739663,10.00812,10.00812005364671] 


Jim FitzSimons 
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December 30 

How do I find the 7 digit number in the range of 10 to 20 that has a base- 
10 log that is closest to the midpoint between two 7 digit numbers? 

That is, I want the log (x, 10) to be of the form: 

1 . abode f_4 999 . . . 
or 

1 .ghijkl_5000 . . . 

I do not care what values a, b, c, d, e, f , g, h, i, j , k, or 1 have. 

But, I do want, mod ( 10 A 6 * log (x, 10) ) to be as close as possible to 0.5 

with x being a 7 digit number. 

For example, 

log (10 . 00706, 10) is 1 . 00030 6_5037 .. . 
log(10. 00918, 10) is 1 . 000398_4 994 . . . 

Some more examples : 
log (10 . 00418, 10) is 1 . 000181_4 971 . . . 
log(10. 00653, 10) is 1 . 000283_5017 . . . 

Fred J. Tydeman Tydeman Consulting 

December 30 

al :=LN (2) /LN (10) * RANDOM ( 1 ) +1 

al : = 1 . 172936545387103 

a2 : = (FLOOR (al*10 A 6) +0 . 5) /10 A 6 

a2 : =1 . 1729365 

a3 :=ROUND (10 A (a2+5) ) /10 A 5 

a3 : =14 . 8 9143 

LOG (VECTOR (a3+x* (0. 00001), x, -1,1), 10) 

LOG ( [14. 8 91 42, 14. 8 91 43, 14. 8 91 44], 10) 

[1.172936112712562,1. 172936404353209, 1.172936695993659] 

Jim FitzSimons 

January 3, 2007 

From Johann Wiesenbauer [ j . wiesenbauer@TUWIEN . AC . AT ] 

Hi Fred, 

Looks like 

APPROX (LOG (15. 51499, 10), 20) = 1.1907515000009954495 

is the best solution to your problem for any x in [10,20] . At least that's 
what my attached DERIVE-program says . 

Cheers , 

Johann 

fred(x := 10, xl := 20, p := 20, r_ := 0.01, s_, t_ := []) := 

Loop 

APPROX (WRITE (x)) 

s_ :: APPROX (ABS (MOD (10 A 6. LOG (x, 10)) - 0.5), p) 

If s_ < r_ 

#1 : Prog 

r_ := s_ 

t_ := AD JOIN(APPROX( [STRING (x) , STRING (LOG (x, 10))], p) , t_) 
x :+ 1/10 A 5 
If x > xl 

RETURN ADJOIN(["x", "LOG(x, 10)"], REVERSE (t_)) 


* 2 : 


fred() 
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*3: 


X 

LOGCx, 10) 

10.0003 3 

1.0000165023767596509 

10.00413 

1.0001314971631705636 

10.006 

1.0002604935473903463 

10.00913 

1.0003934994511117030 

10.01439 

1.0006245005390073631 

10.02999 

1.0013005000247063514 

10.03632 

1.0015744999939723413 

10.17615 

1.0075335000036467535 

10.49053 

1 . 0207994999937194000 

. 15.51499 

1.1907515000009954495 . 


#4: APPROX (LOG (15. 5 1499, 10), 20) = 


1.190751500 


January 3, 2007 

Thank you. That is what I was looking for. 

I see that your program did a brute force try all numbers. 

Now, to see if I can come up with some faster search method. 

Fred J. Tydeman Tydeman Consulting 

Just before finishing this DNL an interesting question came in how to representing the “Klein Bottle” 
together with normals to the surface: 


. . . This is quite a complicated modelling task as not only do I need to de- 
termine the normals to the surface at a given z value, but also to make the 
normal vectors a small enough length to be instructive. Hopefully someone 
can provide some guidance on this issue as I would be looking to apply the 
technique to other surfaces. 



I'll come back to the two forms of the Klein Bottle in the next DNL, Josef 



p 


6 


Bjorn Felsager: Challenges from Fermat 


D-N-L#64 


Solving some Challenges from Fermat 
Using Symbolic Algebra 

See how to use Newton’s geometrical methods to deal 
with some famous number theoretical problems investigated by Fermat. 

Bj 0 rn Felsager; Haslev Gymnasium & HF; Denmark 

Fermat and Diophantine equations 

Back in the 17’th century Fermat acquired an edition of Diophantus, a classical text on number-theory. 
While reading it he made a lot of comments in the margin, notes that were finally published by his son 
in another famous edition of Diophantus. Apart from these notes, Fermat’s work with number theory 
can be followed through his correspondence with many contemporary mathematicians. But these let- 
ters seldom contained many details. In fact only one proof has survived from Fermat - at the very end 
of his copy of Diophantus he found sufficient room to disclose at least some of the details. 

This does not mean that Fermat did not pay attention to proofs. In a famous critique from 1657 of a 
contemporary English mathematician Wallis he comments upon Wallis use of ordinary induction, i.e. 
establishing a hypothesis by gathering a lot of examples supporting a definite pattern, as opposed to 
mathematical induction: 

' One might use this method if the proof of some proposition were deeply concealed 
and if before looking for it, one wished first to convince one self more or less of its 
truth. But one should place only limited confidence in it and apply proper caution. 

Indeed, one could propose such a statement, and seek to verify it in such a way, that it 
would be valid in several special cases, but none the less false and not universally 
true, so that one has to be most circumspect in using it. No doubt it can still be of 
value if applied prudently, but cannot serve to lay the foundations for some branch of 
science, as Mr. Wallis seeks to do, since for such a purpose nothing short of a demon- 
stration is admissible. ’ 

Similarly Fermat writes to Cherselier in 1662: 

The essence of a proof is that is compels belief 

So clearly Fermat was aware of the importance of proofs. This did not prevent him from making mis- 
takes, some obvious, some subtler. But by and large, if Fermat corresponded about something he 
claimed he knew how to prove, it should be taken seriously - despite the fact, that only one proof was 
actually worked out in writing in sufficient detail, to make a reconstruction of his reasoning possible. 

Fermat is most famous for his negative propositions, i.e. theorems that certain problems have no solu- 
tions. E.g. he challenges Sainte/Croix in a letter from 1638 to find 
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• A right-angled triangle, whose area is a square. 

• Two cubes whose sum is a cube. 

• Two fourth powers whose sum is a fourth power. 

And he repeats these challenges to Frenicle in 1640. These problems were solved using the method of 
infinite descent, which he described very briefly in the following way: 

’There is no right angled triangle in numbers whose area is a square ...if the area of 
such a triangle were a square, then there would be a smaller one with the same prop- 
erty, and so on, which is impossible. ’ 

The idea is very ingenious, but without details it is certainly not easy to fill in the gaps in the argu- 
ment! 

But Fermat was also involved in affirmative problems , which he actually found harder. In solving such 
problems he used the method of ascent , which was an extension of a method already described in Dio- 
phantus and known as ’Methodus vulgaris ’. He was very fond of this method: 

which has astonished the greatest of experts, in particular Mr. Frenicle ’. 

In fact he illustrated his method using two examples he learned from Frenicle, who sent them to him in 
1641: 

1) Find a right-angled triangle such that a > b and such that (a - bf - 2b 2 , is a square. 

2) Find a right-angled triangle such that both c and a + b are squares. 

We shall comment on these problems in considerable detail later on, but for the moment we just make 
the remark, that at first Fermat was unable to do anything about them. As he wrote to Carvaci: 

7 saw no way of even tackling them ’. 

But a few years later he can solve them both. In 1643 he communicated the latter solution to 
Mersenne, the triangle being given by: 

a = 1061652293520 b = 456486027761 c = 4687298610289 

In 1644 he reports the solution of the first problem to Frenicle 

a = 1517, b= 156 and c = 1525 

So the years from 1641 to 1643 were very productive for Fermat, and that’s where he laid the founda- 
tion for his favorite and impressive method of ascent. It is this method that we will try to explain in the 
following, but apart from the fact that Fermat revealed no details, it seemed that Fermat solved them 
using purely algebraic techniques, which can be somewhat difficult to follow. Especially since Fermat 
did not possess an efficient algebraic notation and therefore had to spell his arguments in long chains 
of words. In stead we will follow an idea from Newton, and show how to solve such problems using 
geometrical insight. 





Bjorn Felsager: Challenges from Fermat 


D-N-L#64 


P 8 


Newton’s insight 

Newton acquired an edition of Diophantus in 1670, which may in fact be the very one published by 
Fermat’s son thus containing the marginal notes of Fermat. Newton wrote his own comments upon the 
method of Diophantus, and especially he showed how to interpret these methods in geometrical terms. 
He established thereby that: 

• Many number theoretical problems can be converted to a geometrical form dealing with ra- 
tional points on simple algebraic curves. 

• The simplest such curves, i.e. the conics , possess rational parameterizations, which can be 
used to summarize all the rational points. 

• The next simplest curves, the so called elliptic curves , had rules for generating new rational 
points from old ones, which were summarized in the so-called chord-tangent- method. 

These methods we will try to throw some light upon in what follows. 

Pythagorean triples 

Diophantine equations have a long story. One of the oldest problems in number theory is no doubt the 
problem about Pythagorean triples , i.e. triples of whole numbers (a, b , c ) which fulfill the Pythago- 
rean relation a 2 + b 2 = c 2 , i.e. a, b and c are the legs and the hypotenuse of a right angled triangle. 
Specific solutions like the famous (3, 4, 5) has apparently existed back to the earliest records of 
mathematical activities. General rules for findings such triples go back at least to Pythagoras and 
Plato, but they are presumably much older. E.g. the Pythagoreans knew that the sum of the first n odd 
numbers is a square. In modern notation we therefore have: 

1 + 3 + 5 + . . . + (2/7—1) = n 2 . 

But it then follows that two subsequent square numbers have to satisfy the relation: 

n 2 + (2«+l) = («+l) 2 

Provided we can come up with an odd number 2n+l which is also a square t 2 we therefore generate a 
triple, since we now get: 

n 2 + t 2 = («+ 1) 2 , with n = l A(t 2 - 1) . 

Notice, that t is necessarily an odd number. 

When t = 3 we get n = 4, i.e. we generate the triple 3 2 + 4 2 = 5 2 . 

When t = 5 we get n = 12, i.e. we generate the triple 5 2 + 12 2 = 13 2 . 

When t = 7 we get n = 24, i.e. we generate the triple 7 2 + 24 2 = 25 2 . 

But there are other triples, where the hypotenuse does not differ by one unit from the longest leg, e.g. 
8 2 + 15 2 = 17 2 . So how do we generate all Pythagoreans triples? That’s the problem we will now 
solve using Newton’s method! 
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Rational points on the unit circle 

The first observation is that we can rearrange the Pythagorean identity in the following way: 


a 2 + b 2 = c 2 


=> 



= 1 . 


Pythagorean triples ( a , b , c ) therefore generate rational points (x, y) = ( a/c , b/c) on the unit circle : 

x 2 +y 2 =l . 

But this also goes the other way, since we can always find a common denominator for the rational 
coordinates of a rational point. I.e. we can always write it in the form (a/c, b/c), with a, b and c conse- 
quently constituting a Pythagorean triple. This is a very important observation: Typically a Dio- 
phantine problem dealing with integers can in this way be turned into a problem dealing with rational 
points on a simple algebraic curve. 

The next observation deals with rational points on a quadratic curve, i.e. a conic section. It follows 
the following principle, known as the principle of all or nothing : 


Either there are no rational points or they are all given by a rational parameterization. 

We will show this in detail for the unit circle, and then comment upon the general principle after- 
wards. 



The unit circle contains four trivial rational points, the four corner points on the axes: (±1, 0) and (0, 
±1). We will select one of them as the unit point E for the unit circle. Traditionally the corner point E 
= (1,0) is selected. Now consider the straight lines passing through E. They give rise to the following 
nice observation: 


A secant passing through E will intersect the unit circle in another rational point pre- 
cisely when the slope of the secant is rational. 


Proof. First the easy direction. Assuming the secant will intersect in the additional rational point 
(xo, To) the slope will be given by 


*0 - 1 


which will therefore also be rational. 
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2 2 2 2 

This corresponds to the Pythagorean triple: ( u — v , 2 uv, u +v ) . 

Remark: Notice that Pythagorean triples that have a common factor, i.e. are reducible will correspond 
to the same rational points. So essentially we can only hope to produce primitive Pythagorean triples, 
i.e. irreducible triples without a common factor. Furthermore the generators u and v - besides from 
being relative prime -have to have opposite parity. If both are odd they will generate a triple, where all 
numbers are even and hence reducible. We summarize this observation in the following famous struc- 
ture formula for the primitive Pythagorean triples: 


A triple (< a , b , c ) is a primitive Pythagorean triple precisely when it can be written in 
the form 

(a,b,c) = (u 2 -v 2 , 2 uv, u 2 +v 2 ) 

with the generators (w, v) being relative prime, with u > v and the generators having 
opposite parity. 


The Group structure for the unit circle 

The existence of a rational parameterization is closely associated with the group structure of the ra- 
tional points on the unit circle. On the unit circle you can add angles and hence you can add points. 
The interesting point is that you can define this adding procedure in a purely algebraic fashion , that 
furthermore can be generalized to other families of simple algebraic curves, the so called elliptic 
curves. Notice first that the secant connecting two points P u and P v has to be parallel with the secant 
connecting the unit point E and the sum of the points P u+V . In particular they have to have the same 
slope: 



Hence we can define the sum of the two points P and Q as the intersection point between the unit 
circle and the secant passing through the unit point E , having the same slope as the secant PQ. 
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This is where the rational parameterization comes into the game. The sum P + Q will have the slope 

y? — Ti 

relative to the unit point E , i.e. the parameter t from the rational parameterization has pre- 

x 2 - x, 

cisely the rational value: 

t=y^n. 

X 2 ~ X, 

So the group structure on the unit circle is respecting the rational points. 

Substituting the above parameter into the equations for the parameterizations actually makes it possi- 
ble to compute the coordinates for the point P + Q in a purely algebraic fashion, although CAS- 
programs will have great difficulties reducing the end result to its most simple form. 



We will therefore restrict ourselves to a discussion of the duplication formula , where we simply dou- 
ble the angle for a point P(x 0 ,yo) on the unit circle. The tangent is perpendicular to the radius, so the 
slope of the tangent is given by 

x 

ci = 

tan gent 


But since the slope of the tangent at the point P = (x 0 , yo) is given by - x 0 /yo the duplicate point 2 P 
must correspond to the parameter value t = -xo/yo . Consequently the duplicate point IP is given by: 


(x,y) 


' t 2 - 1 -It ' 


( 2 'S 

f-1 - 2 ■ (-—) 
y 0 y 0 


' x o 2 -y 0 2 2x 0 y 0 

[t 2 +iY + iJ 


2 ’ 2 

X \ +1 \ + l 

Ao Jo J 


2 2 ’ 2 2 
Yo + Jo x o +Jo J 


But since the point P is a point on the unit circle we know that x 0 2 + yo reduces to 1 . Hence the for- 
mula for 2 P simplifies to: 


2 P = (x, y) = (x 0 2 -y 0 2 , 2x 0 y 0 ) 
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Once again we can throw in a rational point (x, y ) = (ale, b/c ) corresponding to a primitive Pythago- 
rean triple (a, b, c ). Consequently the duplicate point 2 P gets the coordinates: 

( a 2 -b 2 lab \ 


This corresponds to the primitive Pythagorean triple ( a 2 - b 2 , lab , c 2 ) and thus we have also 
found a duplication formula for primitive Pythagorean triples! E.g. the primitive Pythagorean triple (3, 
4, 5) is duplicated to (4 2 -3 2 , 2*3*4, 5 2 )= (7, 24, 25). Notice the similarity with the structure for- 
mula: The generators for the duplicate triple are precisely the legs of the original triple. It follows that 
you can reverse the process precisely when the hypotenuse is a square. 

This observation makes primitive Pythagorean triples with a square hypotenuse so much more inter- 
esting! 


Fermat’s cubic problem 

Let us now turn to Fermat’s cubic problem dealing with splitting two cubes into a cube, i.e. this time 
we are looking for integers a , b and c fulfilling the relation 

a 3 + b 3 = c 3 . 


We will try solving this problem with a similar technique. The first step is a rearrangement of the 
equation to the form 


a 3 +b 3 = c 3 


M 

3 

+ 

\ b ) 

UJ 


UJ 


= 1 . 


Non-trivial solutions (« a , b , c ) to the cubic Fermat equation thus corresponds to non-trivial rational 
points (x, y) = ( a/c , b/c) on the cubic curve 


* 3 +y 3 = l • 

If it were possible to find a rational parameterization of this curve we would of course control all the 
rational solutions. We could therefore try the same idea as before. There are two trivial rational points 
on the Fermat curve: (1,0) and (0,1). 
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Looking for the secant passing through (1,0) we could therefore investigate additional intersection 
points. But unfortunately there will in general be two additional intersection points. And even if the 
slope is rational - forcing the resulting cubic equation to have rational coefficients - that only allows 
us to conclude, that the sum of the three roots is rational - the two additional roots themselves need 
not be rational. To guarantee an additional rational intersection point, it is therefore necessary to know 
two rational points on the cubic curve in order to compute the additional intersection point with the 
secant passing through these two points. 

But we already have two rational points (1,0) and (0,1). So why not combine these to another rational 
point? Unfortunately it does not work. The equation for the secant becomes x + y = 1. Substituting 
this into the equation for the Fermat curve the cubic equation unfortunately degenerates into a quad- 
ratic equation: 


X 3 +/ = l^x 3 +(l-x) 3 = l^-3x + 3x 2 =0 


(corresponding simply to the two unit points , i.e. x = 0 i.e. x = 1). 

We could also be clever and use one of the unit points as a double point on the curve, i.e using the 
tangent passing through the point. That would also have an equation with rational coefficients and 
hence intersect in a third point with rational coordinates. In principle we need therefore only know one 
rational point, because we can the shoot along the tangent to produce another rational point! 

But once again the idea fails. Both unit points are actually inflection points. As a consequence they 
correspond to triple points on their tangents, so they do not produce additional solutions. What worked 
with such a great success for the unit circle completely fails for the cubic Fermat curve! 

We can however easily find a duplication formula for the cubic Fermat curve. It is done in exactly the 
same way as for the unit circle. This time the duplicate point 2 P is defined as the additional intersec- 
tion point with the tangent through P. It will make it a little easier to symmetrise the problem. So in- 
stead of looking at the equation 


a 3 + b 3 = c 3 . 

we will move the last term to the same side as the other and change the sign of c. The problem then is 
to find integers a , b and c , such that the sum of their cubes vanishes: 

a 3 + 6 3 + c 3 = 0 . 

The corresponding rational point then gets the equation 

x 3 + y 3 + 1 — 0 . 

In this equation we can isolate y to get the equation on its functional form 


j = Ux 3 - 1 
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We are now going to investigate the particular problem of finding a Pythagorean triple where both the 
hypotenuse and the sum of the legs is a square\ First we must reformulate the problem into a geomet- 
ric setting. We start by observing, that we can write both c and a + b as squares, whereas we know 
nothing about the difference a -b: 

c=p 2 , a + b = q 2 and a - b = r. 

So the Pythagorean triple (< a , b , c ) has been replaced by the triple (p, q , r) which satisfies the extended 
Pythagorean identity : 


(a + b) z +(a-bY =2a 2 +2b 2 =2c 2 => q 4 +r z =2p 


r 2 =2p 4 -q 4 


Dividing by q this can be rearranged in the usual manner: 


= 2 — -1 


j = 2x 4 - 1 with j = — and x = — 

q q 


However from the triangular inequality we also get the following constraints 

a-b<c<a + b <=> r<p 2 <q 2 . 

From these it follows that both x and y have to be less than 1 : 0 < x, y < 1 . 

If on the other hand we can find a rational point (x, y) on the algebraic curve y 2 = 2x 4 -1 satisfying 
the constraints 0 < x, y < 1 then we can put (by a suitable choice of denominators) 

p r 

x = — and y = — . 

q q 

Hence the corresponding Pythagorean triple (< a , b , c ) given by 

2 2 

p 2 =c , q 2 = a + b and r = a-b i.e. c = p 2 , a = and b = . 

2 2 

This is the key to finding a solution to the challenge proposed by Frenicle! 


If (< a , b , c ) is a primitive Pythagorean triple satisfying the conditions: 

• The hypotenuse c is a square. 

• The sum of the legs, a + b, is a square. 

Then we can find a rational point (x, y) on the elliptic curve y 2 = 2x 4 - 1 satisfying the 

constraint 0 < x,y < 1 . 

The correspondence is given by 

2 q 2 +r q 2 -r p r 

c = p , a = — - — ,b = — - — and x = —,y = — 

2 2 q q 
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Remark : If instead you tried to find a Pythagorean triple (< a , b , c) where both the hypotenuse c and the 
difference between the legs a - b are squares you can similarly show that they correspond to a rational 
point (x, y) on the same algebraic curve y 2 = 2x 4 -1 . But this time it satisfies the constraint x, y 
>1 . 

Remark: To see why the latter problem is equivalent with Frenicle’s first challenge, we notice, that 
because the hypotenuse is a square the generators (w, v) for the Pythagorean triple {a, b , c ) actually 
form the legs of another right angled triangle ( u , v,p) obtained by bisection!): 

2 2 2 2 2 

u -v =a , 2 uv = b, u +v =c = p 

But the difference between the original legs should also be a square, i.e. the generators must satisfy the 
relation: 

q 2 = a-b = u 2 -v 2 -2 uv = (u-v) 2 -2v 2 
So the root triangle ( u , v,p) precisely solves Frenicle’s problem! 

A summary of Newton’s chord-tangent method 

As we have seen several Diophantine problems end up studying rational points on simple algebraic 
curves. It pays to classify such curves in the following way: 

If the algebraic curve has an equation of the form 

y 2 =p\(x) (where p x denotes a polynomial of degree 1) 
or 

y 2 =p 2 (x) (where p 2 denotes a polynomial of degree 2) 

it is a conic section (in fact it is a parabola in the first case and an ellipse or a hyperbola in 
the second case). 


7 

\ 

y 

y 2 =jjx+l 


X 


1 1 S 




7 

\ 

y 

y 2 =-x 2 +x+4 

\ x 

V 






If the algebraic curve has an equation of the form 

y 2 =p 3 (x) (where p 3 denotes a polynomial of degree 3) 
or 

y 2 = pfx) (where p 4 denotes a polynomial of degree 4) 
it is called an elliptic curve (cubic if the degree is 3 and quartic if the degree is 4). 
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Remark : The name is a little peculiar, since the elliptic curves are not directly related to ellipses. It 
turned out however that curves from the above family naturally pops up if you investigate the arc 
length of an ellipse - hence the name! 

Diophantine problems related to rational points on a conic section are considered trivial. There is a 
complete theory for how to solve such problems: Either there are no rational points or the conic sec- 
tion has a rational parameterization comprising all the rational points. It all follows along the same line 
as we have been discussing in the case of the unit circle. 

Diophantine problems related to rational points on an elliptic curve are considered hard - but within 
reach. There exists a huge theory for the behavior of such points although it is by no means complete. 
Provided the elliptic curve is not degenerate (i.e. does not contain cusps or double points or isolated 
points) it will not have a rational parameterization. But there do exist duplication formulas as well as 
addition formulas which will allow us to systematically build up all the rational solutions. It is pre- 
cisely this structure that is revealed in the chord tangent method going back to Newton. 

Diophantine problems related to rational points on algebraic curves of even higher degree (or rather 
higher genus) are considered extremely difficult. Almost nothing is known of their general behavior. 


Cubic elliptic curves 

We will start by examining cubic elliptic curves y 2 = ax 3 + bx + ex + d . We can then try to combine 
rational points on the elliptic curve thus producing new rational points. This requires two rational 
points to begin with: P\(xuy\) and Pi&i-yi)- Looking at the chord passing through these two points y 
= Ax + B , we notice that it gets a rational slope 


The usual formula for the intercept B 


a = ^A 

x 2 -x 1 

B =y\ — ax i 


shows that B is rational as well. So the equation for the chord has rational coefficients. We now inter- 
sect the chord with the elliptic curve y 2 = ax 3 + bx 2 + ex + d . This will lead us to a cubic equation 

(Ax + B ) 2 = ax 3 + bx 2 + ex + d 
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and since we already know two solutions x\ and x 2 there has to be and additional root of this equation 
x 3 . But the sum of the roots is clearly rational. Hence the third root is also rational. Substituting x 3 into 
the chord equation we find the y-coordinate y 3 = Ax 3 + B . Hence y 3 is also rational. So using the 
chord we have succeeded in combining the two rational points P\ and P 2 into another rational point 
P 3 =Pi& P 2 . 


Should the two rational points be identical i.e. represent a double point all we have to do is to replace 
the chord with the tangent passing through this point. 


An example from Diophantus 

We will exemplify Newton’s method using the following problem from Diophantus (but stating it in a 
modern language!) 


Decompose 6 as a sum of two rational numbers, such that their product is a cube minus its root. 


Decomposing 6 as the sum of 3 + y and 3 - y , we must therefore be able to write their product as 
x 3 -x, i.e.: 


(3 + y ) (3 -y) = x 3 -x =^> 9-y 2 = x 3 -x =^> y 2 = -x 3 +x + 9 

We are thus looking for non-trivial rational points satisfying the equation y 2 = -x 3 + x + 9 as well as 
the constraints x > 1 and |y| < 3 . 



Notice the trivial solutions: (±1,±3) and (0,±3). None of these are admissible by themselves, but they 
can be used to generate admissible solutions. For this purpose we have to use chords or tangents. In 
this case chords will not work, so we will use tangents to generate solutions satisfying the appropriate 
constraint. Looking at the tangent associated with the rational point (-1,3) we succeed: 
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Quartic elliptic curves 

Finally we turn to the case of quartic elliptic curves 

y 2 = ax 4 + bx 3 + ex + dx + e . 

This time a chord is of no use since we will need three rational points if we want to generate an addi- 
tional rational point. In stead we use the parabola passing through the three rational points 

Pi(x u yi\ P 2 (x 2 ,y 2 ) andP 3 (x 3 ,y 3 ). 

As is well known there will be precisely one parabola y = Ax + Bx + C passing through three points 
provided these points do not fall on a straight line. Ignoring these exceptions we proceed with the pa- 
rabola. The coefficients A, B and C can be found solving the following linear system of equations 

yi = Ax 2 + Bx i + C 
y 2 — Ax 2 + Bx 2 + C 
y 2 — Ax 2 + 5x 3 + C 
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Hence the coefficients are rational. 

Intersecting this parabola with the elliptic curve y 2 = ax 4 + bx 3 + cx + dx + e we this time obtain a 
quartic equation (Ax 2 + Bx + Q 2 = ax 4 + &x 3 + cx 2 + dx + e . Since we already knows three solutions 
Xi, x 2 and x 3 to this equation, there has to be a fourth solution x 4 . And since the sum of the roots is ra- 
tional, the same must be true for the last root. Substituting x 4 into the equation for the parabola we get 
the corresponding y-coordinate: y 4 = Ax 4 2 + Bx 4 + C . Hence y 4 is also rational. Proceeding in this 
way we have succeeded in combining the three rational points P u P 2 and P 3 into a fourth rational 
point P 4 = P\& P 2 & P 3 . 

Should some of the given points be identical we just need to interpret the parabola suitably. With a 
double point , the parabola will have the tangent in common with the elliptic curve, i.e. they will have 
the same first derivative. With a triple point the parabola will further more have the curvature in 
common with the elliptic curve, i.e. they will have the same second derivative. 

So if we only know a single rational point P 0 (*o ? yo ) we use it as a triple point. Provided the elliptic 
curve has the tangent slope yd (necessarily being rational) and the second derivative yd f (also neces- 
sarily being rational) at this point P 0 we thus use the parabola with the equation: 

y = '/zjo" ( x - *0 ) 2 + Jo' (x - Xq) + Jo 

Intersecting this with the elliptic curve we thus find the triplication point P 0 & Po & Po- In this way we 
can build up a chain of solutions - unless we have bad luck and just keep generating the same points 
over and over. 


Solving Fermat’s challenge 

As we have seen this deals with the quartic elliptic curve y 2 = 2x 4 - 1 , having the trivial solutions 
(±1, ±1). Due to the high degree of symmetry any solution is actually related to three other solutions 
through reflections in the x- and/or y-axis. This is because the equation is even in both coordinates x 
and y, so that we are free to change the signs. The curve furthermore has the double parabola y = 

± 2 x 2 as an asymptote: 
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But we are really only interested in rational points within the unit square\ However it will be neces- 
sary to combine our way throughout the whole curve. We have four related starting points, so nothing 
will come out of combining them with each other. Instead we will use one of them as a triple point 
P(l,l). So we have to calculate the Taylor polynomial of degree 2, i.e. the quadratic curve passing 
through the same point with the same slope and curvature. But this is easy using a CAS-system. We 
start by isolating;; and then use the Taylor command: 



The result is the quadratic curve 

y = -2(x - l) 2 + 4(x - 1) + 1 = -2x 2 + 8x - 5 . 

We must therefore find the additional intersection point between this parabola and the quartic elliptic 
curve using a solve command. The result is 

x = 13 and y = -239 . 

So the additional point is Q( 13,-239) and we now know four more rational points on the curve. 



Remark. The two illustrations use very different scales, since you cannot put all details into one single 
picture. 

Unfortunately they all fall outside the unit square. We must therefore try again. We will start by 
switching to the point Q\ 13,239) in the first quadrant. We then need to work our way back to the unit 
square. Using Q as a triple point unfortunately does not help. So we will have to combine P and Q. 
This can be done in essentially two different ways: Either we can use P as a double point and Q as a 
single point or vice-versa. It turns out that the first solves the first of Frenicle’s challenges, while the 
second solves his second challenge. This time the parabola therefore only have to be tangential at the 
point Q . We therefore build its equation in the following way 

7 / — 7 i 2 f 8788 ^ __ 57123 

y = taylor(y 2x -1, x, 1, 13) + a*(x-13) = a - x + 26 a -x + 169 a 

^ 239 J 239 
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We then need to adjust the coefficient a so that it also passes through (1,1), i.e. to solve the equation 

, . 2 ( 8788 57123 . . , J 

solve(y = a-x + 26 a -x + 169a ,a)\x = land v = 1 

l 239 ) 239 


24287 

which produces the solution a = . So this time the parabola gets the equation 


_ 24287 -x 2 637 x 8353 

y ~ 17208 8604 17208 


Computing the intersection points between this parabola and the quartic elliptic curve we this time get: 


2165017 

2372159 


and y = 


3503833734241 

5627138321281 


So we did it! We landed right inside the unit square and we have essentially solved Fermats challenge! 
The x- value already gives us the hypotenuse as well as the sum of the legs: 


x=— => c =p 2 = (2 1 650 17) 2 = 46872986 10289 


a + b = q 2 = (2372159) 2 = 5627138321281 

Notice that the denominator of the ^-coordinate precisely coincides with q. From the y-value we then 
obtain the difference between the two legs a and b\ 

a-b = r = 3503833734241 

Hence the legs are given by 


2 2 

a= = 106165 2293520 b= ^ - =4565486027761 

2 2 

The final words 


We will let Fermat have the last word. In a letter to Huygens he concludes: 


‘Such is in brief the tale of my musings on numbers. I have put it down only because I 
fear that I shall never find the leisure to write out and expand properly all these 
proofs and methods. Anyway this will serve as a pointer to men of science for finding 
for themselves what I am not writing out, particularly if Mr. Carvaci and Frenicle 
communicate to them a few proofs by descent that I have sent them on the subject of 
some negative propositions. Maybe posterity will be grateful to me for having shown 
that the ancients did not know everything ' 
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Lower and upper sum: 


At 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

number 

60 

50 

50 

160 

390 

440 

675 

875 

990 

820 

640 

540 

At 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

number 

500 

475 

475 

610 

850 

725 

475 

300 

240 

160 

110 

100 


^/(x,.)-Ax,. =l-(60 + 50 + 50 + 160 + 390 + 440 + 675 + 875 + 990 + 820 + 640 + 

i = 1 

+ 540 + 500 + 475 + 475 + 610 + 850 + 725 + 475 + 300 + 240 + 160 + 110 + 100) = 10235 


At 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Anzahl 

100 

75 

150 

280 

420 

640 

940 

1000 

1020 

990 

825 

640 

At 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Anzahl 

525 

500 

600 

825 

910 

900 

610 

450 

310 

240 

150 

120 


24 

■ Ax,. = 1 • (100 + 75 + 150 + 280 + 420 + 640 + 940 + 1000 + 1020 + 990 + 825 + 

1=1 

+ 640 + 525 + 500 + 600 + 825 + 910 + 900 + 610 + 450 + 310 + 240 + 150 + 120) = 13220 


Average: 


13220 + 10235 
2 


« 11728 


Questions a) - c) didn't need any technology except ruler, paper and pencil. This is now the moment 
when technology comes into play. For the following we read off the “traffic density values” from the 
graph and either put them into a data sheet on the TI or put them into a matrix with DERIVE. I will 
start with the Tl-process: 


Additional question for the students: Right- and Leftsum have the same value? Why can this happen? How to 
obtain these to sums with one single command? 
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d) The scatterplot of the data points shows a form with two peaks and we have the idea that a quartic 
might fit (the TI does not provide polynomial fits of higher degree!). 



It turns out that the quartic regression does not deliver a good model function, but that the value 
for the integral seems to give a useful estimation for the total number of vehicles. 


e) and f) Cannot be worked through with this bad model - see the DERIVE-version and the TI- 
extension. 


g) If we accept the quartic being good enough for integrating the sum, then it makes sense to apply 
numeric methods: Midpoint- and Trapezium Method (left) and Simpson Method (right). 
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We must still be aware of the fact that the model function is not the best one! ! 

j) The average number per hour is approximately 490 vehicles. 

k) What do you think? 

l) See the DERIVE procedure. 

But first Ed like to add a possible extension using the Tl-program polreg from DNL#63, which over- 
comes the quartreg-restriction: 



The polynomial regression function of degree 6 looks a bit better, but the estimation of the number of 
vehicles (area under the curve) is obviously worse. We try a polynomial of degree 8: 



This is a again bit better, but it must be said that we did not have polreg at our disposal when we de- 
signed the example two years ago ... 

Suddenly we had the idea that each of the peaks reminds us on the bell curve of the normal distribu- 
tion. We took again the TI and tried to find a bell for the morning and another one for the afternoon by 
varying the parameters ... 

In most of the statistics textbooks you can find a sketch of a “bimodal distributions” among the con- 
tinuous distributions, but without presenting any example! Here is one. 


I did the same with DERIVE varying the parameters using the slider bars. This is a short report on the 
DERIVE procedure: You see the last two rows of matrix veh (the data points from above) and the 
polynomial fit of degree 10: 

23 120 

. 24 100 . 

Answers a) - c) are identical to the Tl-version 
We try a polynomial of degree 10 fit: 

[ 10 98765432 1 

#4: pollO := FIT(Lx, a»x + b»x + c»x + d*x + e»x + f»x + g»x + h»x + i»x + j*x + kj , veh) 

-7 10 -5 9 3 7 6 

#5: pollO := 4.3431751.10 »x -6.0325023.10 .x + 0.0031706072.x -0.091235174.x + 1.5637329.x - 

5 4 3 2 

16.213440.x + 99.334364.x - 339.55054.x + 600.33404.x - 423.20810.x + 109.25742 

#6: TABLE(poll0, x, 0, 24, 0.01) 

d) The interesting part begins now. First of all see some polynomial regression curves: 
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Now let's follow the “Gau 6 Curve Idea”: I define the generalized “Morning Traffic Bell” and intro- 
duce slider bars for the parameters a , b , 5 and d. It is a fine question for the students to estimate the 
values for the parameters (vertical stretch factor a , the mean b , the standard deviation s - they might 
know that the distance between turning point = mean and inflection point is the standard deviation 5 - 
and d is some vertical shift). Having this in mind it should not be too difficult to install meaningful 
ranges for the parameters. 

2 

- 1/2. ((x - b)/s) 

m_tr(t) := a + c + d 



This is the Morning Traffic Curve. In a similar way we can find an Aftemoon/Evening Traffic Curve. 
Then we have two bell curves. How to obtain a curve for the whole day? Simply add the functions. As 
you can see the afternoon part of the Morning Bell has little influence on the sum curve and the morn- 
ing part of the Evening Bell has only very little influence on the sum curve, too. So we get the vehi- 
cle_density curve for the whole day. Just play a little with the slider bars to find the best fit and it 
is done. 
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2 


2 


*13: - 1/2. ((x - S. 3)/2.5) - 1/2. ((x - 16. 8)/2. 4) 

veh_dens(t) := 962* c + 87S*c +32 

I was just interested how MS-Excel would solve the problem using the solver to minimize the SSE 
(Sum of Squared Errors). This is the plot created by means of the Excel Solver: 



I compared the sum of the squared errors of „my“ double bell shaped-fit-curve and the Excel-Curve of 
best fit: 

2 2 


#14: - ((t - 8.545)/3.526) - ((t - 16. 674)/3 . 025) 

excel_dens(t) :i 90S.16.e + 726. S*e + 88.33 


*15: I VECTOR 


ijv 

5: l|v 


)) 


*16: I VECTOR 


(veh - veh_dens(veh )) , k, 1, DIM(veh) 
k, 2 k, 1 


(veh - excel _dens (veh )) , k, 1, DIM (veh) 

k, 2 k , 1 


= 135964.24 


60312.69 


And I found out that - thanks Bill Gates - that this best fit curve was much better - even without 
reaching the peak in the afternoon! 






p30 


Thomas Himmelbauer: A Traffic Density Problem 


D-N-L#64 


I find it very challenging and useful to switch between the tools - even in school. Most students are 
accustomed to a spread sheet and why not change to it and use its results in CAS and sometimes the 
other way round? The students should learn to become flexible in using their tools, including brain, 
paper and pencil! 


It is not forbidden to become cleverer, so we will use the Excel-fit-curve for further investigations. 



We will do the further investigation using the best fit and provide a shorter name for it, trd(t) 

2 2 

#19: - (ft - 3 . 545)/3 . 526) - ((t - 16. 674)/3 . 025) 

trd(t) := 90S.16*e + 726. S*c + S3. 33 


# 20 : 

# 21 : 

# 22 : 


N50LVE 


— trd(t) = 0, tj, N50LVe|— trd(t) = 0, t, 10, uj, N50LVe|— trd(t) = 0, t, 14, 2ojJ 


[t = 8.5515292, t = 12.938050, t = 16.635525] 



' 3.55 

997.02402 ' 

TABLE(trd(t), t, [3.55, 13, 16.64]) = 

13 

433.60933 


. 16.64 

319.70609 . 


The local maximum values are reached at 8.6 in the morning and 16.6 in the afternoon with a 
density of approx 1000 and 820 vehicles per hour. Minimum traffic is at 1 pm with about 440 ve- 
hicle per hour. The nightly minimum cannot be calculated as a local extremal value. According to 
the data points it is approximately at 2 am, the model offers a boundary extremal value a 0 
o’clock. 


f) We find and interpret the values of the 1 st derivative for t = 7 and t = 21 


#23: TABLE (trd r (t) , t, [7, 21]) = 


7 136.33365 

21 -33.904376 


At 7 o’clock traffic density increases by 186 cars/per hour and at 21 o’clock it decreases by 89 
cars/hour. 

g) Total number of all cars passing this day is the area under the curve, i.e. the defined integral. We 
start using some numerical methods in the same way as with the TIs. We can assume that the re- 
spective functions for numerical integration are ready made and the students can load their tool 
box. 
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1 b - a n 

#24: trapsum(u, v, a, b, n) := I 

2 n i=l 


SUBST 


(■i - l).(b - a) 


u, v, a + 


+ SUBST 


^u, V, 


i . (b - a) 


)] 


#25: midptsum(u, v, a, b, n) := 


b - a n 


#26: simps(u, v, a, 



n 

#27: 

trapsum(trd(t) , t, 

#23: 

midptsum(trd(t) , t 

#29: 

simps(trd(t) , t, 0 


n f 

I SUBST u, v, 

i=l l 

n b - a f f 

b, n) := I SUBST u, v, a 

i=l 6*n l l 

3 1 f 

+ SUBST u, v, a + 

n J l 


b - a (i - l)*(b - a) 
a + + 


n n 

(i - l)*(b - a) 


n 


+ 4* SUBST 


ju, V, 


(i - 1) ♦ (b - a) b - a 

+ 


i • (b - a) 


]) 


All estimated results are « 11690 vehicles. 


h) Calling finally the defined integral should give a very similar result: 

24 

#30: J trd(t) dt = 1601 . 0860 . Vtt*ERF(4. 3831537) + 1601 . 0860.7*. ERF(2. 4234259) + 1099. 285 •Jn* ERF(5 . 5120661) 

0 


+ 1099. 285. y*.ERF(2. 4218181) + 2119.92 


24 

#31: J trd(t) dt = 11689.557 

0 

Appearance of the error function in #30 might lead to some discussion. 


i) Now we meet the undefined integral. 



k) I hope that the students will find some examples (consumption of water, electricity during a day, 
touristic seasons, . . . .) 


1) The last point is really a challenge and a great occasion to demonstrate the power of a CAS. This 
point is always a final and closing highlight when working through this example with teacher stu- 
dents and with experienced teachers on seminars. 

We want to extend this function over a week - let's say without weekend, because these two days 
will follow other rules. So we have a periodic function and periodic functions provoke the idea of 
involving trigonometric functions. 
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The moment will come when somebody whispers “Fourier??”. That is it, but how to perform a 
Fourier analyse? Most of us - except teachers on technical schools - did this the last time during 
their university studies. 

Easy done let's ask DERIVE. We call the Online Help, type in Fourier into the Index search and, 
here it is: 


FOURIER(y, t, tl, t2, n) simplifies to a Fourier series approximation of the expression y(t) as the 
variable t varies from tl to t2, truncated to the nth harmonic. The Fourier series has the form bO 
plus the sum of 

ak-sin(2-p-k-t/(t2-t1)) + bk-cos(2-p-k-t/(t2-t1)) 

for k=1 to n, where bO through bn and al through an are expressions independent of t. If simpli- 
fying an expression involving the FOURIER function does not provide a closed form, you can try 
approximating it instead. For example, to see a Fourier series approximation for a square wave 
truncated to the 12th harmonic, simplify and plot the expression 

FOURIER (SQUARE_WAVE (t) , t, 0, 2, 12) 

So we use this provided tool: 

The exact result amazes the students, the approximated form look a bit “friendlier”. 

#34: FOURIER(trd(t), t, 0, 24, 4) 

#35: - 9. 2944433 *CQS(1 . 0471975 ■ t) - 13.735045 *5IN(1 . 0471975 *t) + 131 .41071 *CGS(0. 73539316* t) + 

63. 374177 *SIN(0. 73539316* t) - 131 . 12044*CQ5(0. 52359377 *t) - 35 . 037440 *SIN(0 . 52359377 * t) - 
330. 92730*COS(0. 26179933* t) + 39.435113 *SIN(0. 2617993S*t) + 437.06490 



cally. 

One of my colleagues was so astonished that he said: “This is really wonderful, I am feeling like 
an apprentice of a wizard” (“Das ist ja wunderbar, ich fuhT mich wie ein Zauberlehrling”). 
Maybe that students might feel like Harry Potter. 


I will close thanking Thomas Himmelbauer for designing this wonderful example which offers so 
many occasions to work on calculus problems, which forces the students to interpret the results 
and which gives the opportunity to use various tools up to change the technology. He is the au- 
thor of many other examples which can be found an the homepage of the ACDCA. 


www . acdca . ac . at 
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In the last DNL Peter Schofield and I showed some plots of contour- or level curves 
and as a useful side product ways to program implicit 2D-plots. Families of 3D- 
level curves give a nice presentation of surfaces and solids. Peter did a lot more - he 
wrote tools for presenting 3D implicitly given surfaces. I present some results to 
show the use of his functions. Many thanks, Peter. 


Implicit Plots in 3D 

Peter Schofield 


The following three functions do all the work: 


Impl icitPts(f , st := -5, end :i 5, step :i 0.25) 

Impl icitDots(f , st := -5, end := 5, step := 0.25) 
InterPtsff, g, st :i -5, end :i 5, step :i 0.1, rad := 0.3) 


If you are satisfied with the standard settings (-5 <x,y,z< + 5 and a scanning step of 0.25) then you 
only need to enter the function(s) in implicit form(s). For changing the settings you can enter vectors 
for st, end and step. See a demonstration of applying standard settings and “individual settings”: 



Lets start with a simple cone x 2 -y 1 +z 2 = 0 appearing in “Dots” and in “Points” (Size Me- 
dium) 

2 2 2 2 2 2 

Impl icitDots(x - y + z ) Impl icitPts(x - y + z ) 
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The next surface is a product of my fantasy. I compare the 3D-plots created by DERIVE, DPGraph and 
Autograph: 


4 3 2 


Iiip’HcitDotsCx + (y - 1) + (2-z) = 32) 




DERIVE 


DPGraph 



Autograph 

Now we have a torus with special values for scanning ranges in x-, y- and z-direction. 


2 2 2 2 2 

Iinp"HcitDots(36 ■ z + (x + y + z - 10) = 36, [-4, -4, -1], [4, 4, 1]) 
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Peter produced one large sphere surrounded by 8 smaller spheres. For convenience define the 
SPHERE-function and then put them all together: 

2 2 2 2 
sphere(a, b, c, r) := (x - a) + (y - b) + (z - c) - r 


Impl icitDots (sphere (3 , 3, 3, 1), 2, 4) 

Iiipl icitDots (sphere (3 , 3, -3, 1), [2, 2, -4], [4, 4, -2]) 
Impl icitDots (sphere (3 , -3, 3, 1), [2, -4, 2], [4, -2, 4]) 
Impl icitDots (sphere (3 , -3, -3, 1), [2, -4, -4], [4, -2, -2]) 
Impl icitDots (sphere (-3 , 3, 3, 1), [-4, 2, 2], [-2, 4, 4]) 
Iiipl icitDots (sphere (-3 , 3, -3, 1), [-4, 2, -4], [-2, 4, -2]) 
Impl icitDots (sphere (-3 , -3, 3, 1), [-4, -4, 2], [-2, -2, 4]) 


Impl icitDots (sphere (-3 , -3, -3, 1), -4, -2) 
Iiipl Hoi tDots (sphere (0, 0, 0, 2), -2, 2) 



InterPts(surfacel, surface 2, [settings for domain, step and rad]) shows the inter- 
section of implicitly given surfaces. This is Peter's advice: 


Intersections of surfaces 

For each of #36-#39 3D-plot use: Insert > Plot > Fled Hun > Custom 
Select Point Colour using Top (Hi n, Flax) and Dot Colour using Grids. 

Set Apply parameters to rest of plot list (ON) and then Finish (be patient): 
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We produce the intersection of two planes and the intersection of a cone and a plane resulting in a 
hyperbola: 




Z Z Z 

InterPtsfx + Z-y + z = 3, x - y + 7 ■ z = 0) 


InterPtsfx +Z,x +y -z) 

Impl -ici t Dot s(x + Z-y + z = 3) 

#37: 

IinplicitDotsfx + Z) 



Z Z Z 

Inpl-ic-itDots(x - y + 7-z = 0) 


IiiplicitDotsfx + y = z ) 




Sphere and cylinder give in a special configuration the Window ofViviani : 

2 2 2 

Impl icitDots(x + y + z = 16, -4, 4) 

2 2 

Impl icitDots((x - 2) + y =4, -4, 4) 

2 2 2 2 2 

. Inter Pts(x + y + z = 16, (x - 2) + y = 4 , -4 , 4 , 0. 05 , 0.05) . 
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Intersection of cylinders: 



3D_PL0T (notice all three solution planes are plotted): 

2 2 

Impl icitPts((x - 1) *(y + l)*z = 0) 

(If you try this equation with 'DPgraph' you don't get the correct plot?) 


x 



DPGraph 


DERIVE 



Autograph 
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Contour Pts_XY, Contour Pts_XZ, ContourPts_YZ and Contour Dots_XY, Contour Dots_XZ, 
Contou r Dots_YZ were introduced in DNL#63 


Surfaces of Revolution, etc. - try these with Insert>Plot> Large Points > Rainbow! 

Set Apply parameters to rest of plot list (ON) and then Finish: 

2 2 


VECTOR (Con toil r Pt s_XY (x + y 


VECTOR 


r 2 2 

ContourPts_XZ x + z 


9 + 3 -COS(z) , z_, -4, 4, 0.2), z_, -5, 5, 0.2) 


2- COSH 

' y 1 

, y_, -4, 4, O.2I, y_, -5, 5, 0.2 


, 2 j 

J J 



2 2 

VECTOR (Con tour Dot s_XY(x +y - 16 ■ COS (2 -z), z_) , z_, -5, 5, 0.1) 


2 2 2 

VECTOR (ContourDots_XY(x + y + (z - 4) = 16, z_, -4, 4, 0.2), z_, 0, 4, 0.2) 

2 2 

VECTOR (Con ton r Dot s_X Y (x +y = 0.4-(l + 0.1-z), z_, -1, 1, 0.1), z_, -4, 0, 0.2) 

2 2 

VECTOR (Con tour Dot s_XY(x +y = z, [z_, -4], -2, 2, 0.1), z_, 0, 4, 0.8) 

2 2 

VECTOR (Con tour Dot s_XY(x +y = z, [z_, -3.8], -2, 2, 0.1), z_, 0, 4, 0.8) 
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With Peter's special recommendations (one of his favourite graphs): 

Imp! icitPts(SIN(x) + SIN(y) + SIN(z) = 0) 



I found the information in the German edition of the Scientific American (= Spektrum der Wissen- 
schaft) from January 2007 that the Association of German Mathematicians (DMV = Deutsche Mathe- 
matiker-Vereinigung) named a special and spectacular mathematical object after a German writer and 
poet, Hans Magnus Enzensberger (“Der Zahlenteufel”, “The Number Devil”). This is the DERIVE - 
image of this impressive star with six vertices: 


Magnus Enzensberger Star 

222222 222 

Impl icitDots(400* (x *y + y *z + x»z) = l-x - y - z , -1, 1, 0.025, 0.025) 

222222 222 
Impl ic itPts(400 ♦ (x *y + y ♦ z +x ♦z)=l-x -y -z, -1, 1, 0.05, 0.05) 




Magnus-Enzensberger-Star created by DPGraph (setting Color by Steepness) 

- (x + y) 

z = - 4 

2 

VECTOR (ContourDots_XV(LN((z + 4) ) + x + y, z_, -5, 5, 0.5), z_, -3.9, 5, 0.05) 

2 

Imp! ic itPts(LN((z + 4) ) + x + y) 



5 5 5 5 

Impl icitDots( |x | + |y| + |z| - 2 , -2, 2, 0.125) 
5 5 5 5 

Impl icitPts( | x | + |y | + |z| - 2 , -2, 2, 0.1) 
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Morley in the Mirror 

Peter Luke-Rosendahl, Germany 

Given is a triangle ABR. Intersection of reflection of AB wrt AR and reflection of AB wrt 
BR gives intersection point O. Draw two lines which form an angle of 30° with OR and 
intersect them with AO and BO giving points Q and P. 

(1 ) Show that the Morley Triangle PQR is equilateral. 

Further reflections of AR and BR with respect to AQ and BP intersect in point C. It is 
obvious that angles CAB and CBA are divided into three equal parts by the straight 
lines, but 

(2) Show that lines CQ and CP also divide angle ACB into three equal parts. 


Before starting with Derive and the computer algebra realisation and proof I'd like to make clear the 
problem using the free dynamic geometry program GeoGebra. You see the plot and the last part of 
the protocol which is verifying both statements from above. 

■■ 9 

■a 
-a 
■a 
-a 
a 
-a 
■a 
a 
-a 
■a 
a 
-a 
a 

Moreover we can ask for an accurate report on 
the construction process so we can reproduce it 
point for point and line for line. 

You can see how to find point C together with all 
coordinates and equations of lines (or circles or 
conics if there were some to appear). 

I find this as a good starting point to make the 
next step: the formal proof. 

This will be done by Peter's Derive paper. 

Peter worked with fixed points A, B and R and 
did the proof finally with generalized coordinates. 

I follow his tracks but include slider bars for the 
coordinates and have one fixed triple A_, B_, R_ 
for performing parallel numerical checks. 

The screen shots show various basic triangles 
(various settings in the slider bars). 


S Construction Protocol 
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#1: CaseMode := Sensitive 

#2: [bx :=, rx :=, ry :=] 

#3: [A := [0, 0], B := [bx, 0], R := [rx, ry]] 

*4: [A_ := [0, 0], B_ := [IS, 0], R_ := [12, 3]] 


We start by providing a toolbox for further use . 

Line (= Gerade) defined by points U,V with parameter s 

*5: g(U, V, s) : = U + s*(V - U) 

Pedalpoint (= Lotfuftpunkt) of P0 with respect to linefU.V): 

(P0 - U) * (V - u) 

LFP(P0 , U, V) := U + .(V - U) 

* 6 : 2 

(V - U) 


Result of reflecting P0 with respect to line(U,V) (=Bildpunkt): 

#7: BP(P0, U, V) := P0 + 2 • (LFP(P0 , U, V) - P0) 

Intersection point (= Schnittpunkt) of two lines gl and g2: 

#8: SP(gl , g2, p, q) := SUBST(gl, p, FIRST (FIRST (SOLUTIONS (gl = g2, [p, q])))) 

Rotation about a": 


#9: 


rotation (a) :s 


‘ COS (a) 
. SIN(a) 


- SIN(a) ' 
COS(ot) . 


#10: TRIANGLE(U P V, W) := [U, V, W, U] 

#11: TRIANGLE (A, B, R) 

Figure 1 

#12: [BP(B, a, R), g(A, BP(B, A, R)), BP(A, B, R), g(B p BP(A, B, R))] 

Figure 2 




#13: 0(U, V, W) := SP(g(U, BP(V p U, W) p p), g(V p BP(U p V, W) p q)) 

#14: [0(A P B, R) , R] 

#15: gQ(A, B, R, p) := R + p. rotation (3 0° ). (0(A P B p R) - R) 

#16: gP(A p B. R, p) := R + q.rotation(- 30®).(O(A p B, R) - R) 

#17: [gQ(A p B, R, p), gP(A p B, R, p)] 


Figure 3 
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figure 3 figure 4 

Finding the Morley Triangle points Q and P: 


#13: 

P(U, V, W) := SP(gP(U, 

V, 

W> q), g(v, 

> 

o 

> p), q. 

■ P) 

#19: 

Q(U, V, W) := SP(gQ(U, 

V, 

W, p), g(U, 

0(u, V, w) 

. q)> Pp 

■ q) 

#20: 

MorleyTriangle(U, V, 

W) 

:= TRIANGLE (P(U, V, W) , I 

Q(U, V, 

w), w) 

#21: 

norleyTriangle(A, B, 

R) 






Figure 4 

We have to show that all angles of APQR are 60° (#22-#24) and/or that sides PQ, PR and 
QR are of equal length (#25, #26). We measure the angles in positive direction: 


ACOS 


(v 


#22: 





Uv - 

u|. |w - u| ) 


w«(u, 

V, 

W) := 


1° 

#23: 

W«(R, 

QCa, b, 

R). 

P(A, B, 

R)) = 60 

#24: 

W«(Q(A, 

B, R) 

, P(A 

, B, R) 

, R) = 60 

#25: 

|Q(A, 

B, 

R) - 

R| - 

|P(A, 

B, R) - R | = 0 

#26: 

|QCA, 

B, 

R) - 

P(Ap 

B, R) | 

- |P(A, B, R) - R| = 0 


Finding point C: 

Reflection of AR and BR wrt to AQ and AP gives two lines, their intersection point is C 
We calculate and plot lines CQ and CP. 


#27: 

gAC(U, V, W, p) := g(U, 

BP(W, 

u, 0(U, 

V, W))p p) 


#23: 

gBC(U, V, W, q) := g(V, 

BP(W, 

v, o(u, 

V, VO), q) 


#29: 

C(U, V, W) := SPCgAC(U, 

V, W, 

p), gBC(U, V, W, q), 

p, q) 

#30: 

[gAC(A, B, R), gBC(A, 

B, R), 

C(A, B, 

R)] 


#31: 

[g(C(A, B, R), Q(A, B, 

R)), 

g(C(A, B 

1 R)p P(A, B, 

R))] 


Figure 5 
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figure 5 

Finally we want to show that gCQ and gCP form a third of angle ACB with sides AC and BC. 
This is no problem using basic triangle A_B_R_ (#33, #34), but does not work in the general 
form (expression #32 simplifies to a huge expression without any useful meaning). 


#32: IF(3.W«(A, C(A, B, R), Q(A, B, R)) = W<x(A, C(A, B, R), B), Ok, false) 

/33: IF(3.W«(A_, C(A_, B_, R_), Q(A_, B_, R_)) = W«(A_, C(A_, B_, R_) , B) , OK, false) 
#34: IF(bx >0, OK, false) 

So we try the other way round: Q1 is mirror image of Q with respect to line PC. Our conjec- 
ture is true if Q1 is element of line BC. 


#35: Ql (A, B, R) := BP(Q(A, B, R) , P(A, B, R) , C(A, B, R)) 
#36: SOLVE (B = g(C(A_, B_, R_) , Ql(A_, B_, R_) , s)) 

235.^3 235 

#37: bx = 13 A s = 

187 374 

#38: SOLVE (B = g(C(A, B, R) , Ql(A, B, R)), s) 


2 2 2 2 2 
(bx - 2 ♦ bx * rx + rx + ry )*(3*rx - ry ) 


#39: (s = ±» a ry = 0) v 


2 * ry * (rx + ry ) * (*/3 * bx - * rx + ry) 


Expression #37 demonstrates that Ql is on BC for the special case and #39 does the same 
for the general case. The same procedure can be done for point P. 

#40: PI (A, B, R) := BP(P(A, B, R) , Q(A, B, R) , C(A, B, R)) 

#41: SOLVE(A = g(C(A, B, R), PI (A, B, R)), s) 

2 2 2 2 2 
(rx + ry )*(3*bx - 6*bx*rx + 3*rx - ry ) 

#42: s = 


2 2 2 
2*ry*(bx - 2*bx*rx + rx + ry )*(i/3*rx + ry) 





D 


-N-L#64 


Peter Liike-Rosendahl: Morley in the Mirror 


p45 


There are some restrictions for s dependent on the coordinates, because of 3a + 3|3 < 180° 
follows a + p < 60°. 

The intersection points of the angle trisecting lines form an equilateral triangle 


What happens if R is moving on a line parallel to the x-axis? 


#43: HorleyFam(t) := VECTOR (Mori eyTr Hang! e(A_, B_, [t, 4]), t, 0.2, 17.8, 0.2) 



What is the locus of points P and Q? The power of CAS helps to give an answer: 

#44: Q(A_, B_, [t, 4]) 

S*73*t*(73*t - 18*73 - 4) 


32 *i/3 • (i/3 *t - 6*i/3 - 4) 4*i/3 

+ t 


#45: 

3 * (t - lS*t - 24* ,/3 + 16) 
#46: P(A_, B_, [t, 4]) 

32 *,/3 • ( 1/3 *t - 12.i/3 + 4) 


#47: 


3 * (t - 18*t - 24*73 + 16) . 
4*73 32*73 *(t - 4*73) 


+ t + 


3 * (t - 18*t - 8* (3 .73 - 2)) 

#48: VECTOR(Q(A_, B_, [t, 4]), t, 0, IS, 0.01) 
#49: VECTOR(P(A_, B_, [t, 4]), t, 0, IS, 0.01) 


3 * (t - 18*t - 8*(3* 73 - 2)) 


+ 8 



-14 -12 -10 -8 -6 -4 -2 


2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 31 
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The locus for both points with -100 < t < 100 has an interesting form: 



A nice Morley family picture at the end 


HorleyFam(t) := VECTOR 


MorleyTriangle(A_, B_, [t, 2*SIN(t)]), t, - 2*71, 2*71, 


71 ' 

40 > 




